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Abstract 

We study phase diagram of the dense holographic gauge matter in the Sakai-Sugimoto model in 
the presence of the magnetic field above the deconfinement temperature. Even above the decon- 
finement, quarks could form colour bound states through the remaining strong interaction if the 
density is large. We demonstrate that in the presence of the magnetic field for a sufficiently large 
baryon density, the multiquark-pion gradient (MQ-yi^) phase is more thermodynamically preferred 
than the chiral-symmetric quark-gluon plasma. The phase diagrams between the holographic mul- 
tiquark and the chiral-symmetric quark-gluon plasma phase are obtained at finite temperature and 
magnetic field. In the mixed MQ-V¥? phase, the pion gradient induced by the external magnetic 
field is found to be a linear response for small and moderate field strengths. Its population ratio 
decreases as the density is raised and thus the multiquarks dominate the phase. Temperature de- 
pendence of the baryon chemical potential, the free energy and the linear pion gradient response of 



the multiquark phase are well approximated by a simple analytic function wl — ^ inherited from 
the metric of the holographic background. 
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I. INTRODUCTION 

Discovery of the AdS/CFT correspondence jl| and the generalization in terms of the 
holographic principle have provided us with alternative theoretical methods to explore the 
physics of strongly coupled gauge matter. Holographic models have been constructed to 
mimic behaviour of the strongly coupled gauge matter in various situations. The Sakai- 
Sugimoto (SS) model (2), is a holographic model which contains chiral fermions in the 
fundamental representation of U(N C ). Its low energy limit is the closest holographic model 
of the QCD so far. It can also accommodate distinctively the chiral symmetry restoration 
and the deconfinement phase transition in the non- antipodal case (4]. It provides interesting 
possibility of the existence of the exotic nuclear phase where quarks and gluons are deconfmed 
but the chiral symmetry is still broken. 

In the SS model, there are two background metrics describing a confined and a deconfmed 
phase. The deconfmed phase corresponds to the background metric with a black hole horizon. 
The Hawking temperature of the black hole is identified with the temperature of the dual 
"QCD" matter. When gluons are deconfmed, the thermodynamical phase of the nuclear 
matter can be categorized into 3 phases, the vacuum phase, the chirally broken phase and 
the chiral-symmetric phase. In the deconfmed phase, the interaction between quarks and 
gluons become the screened Couloub potential. If the coupling is still strong, bound states 
of quarks could form (see Ref. js-lO] for multiquark related studies). The phase diagram 
of the holographic nuclear matter in the SS model is studied in details in Ref. Uj and 
extended to include multiquarks with colour charges in Ref. 10j. It has certain similarity to 
the conventional QCD phase diagram speculated from other approaches e.g. the existence of 
critical temperature line above which chiral symmetry is restored. The phase diagram also 
shows the thermodynamic preference of the multiquark phase with broken chiral symmetry 
for moderate temperature in the situation when the density is sufficiently large. As an 
implication, it is thus highly likely that matters in the core of neutron stars are compressed 
into the multiquark nuclear phase. A thorough investigation on the multiquark star suggests 
higher mass limits of the neutron stars if they have multiquark cores 121 ] . 

When the magnetic field is turned on, the phase structure becomes more complicated. 
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Magnetic field induces the pion gradient or a domain wall as a response of the chiral con- 
densate of the chirally broken phase 13]. In the confined phase, this is pronounced fl^ |. 
However, it is demonstrated in Ref. [ljthat the pion gradient is subdominant to the contri- 



bution 
Ref. 



rom the multiquarks in the chirally broken deconfined phase. It was also shown in 
15| that for sufficiently large density, the multiquark phase is more thermodynamically 
preferred than the chiral-symmetric quark-gluon plasma for small and moderate magnetic 
field strengths. Therefore it is interesting to explore the phase diagram of the deconfined 
nuclear matter in the presence of the external magnetic field. We establish two phase dia- 
grams between the chirally broken multiquark (xSB) and the chiral-symmetric quark-gluon 
plasma (xS-QGP), one at fixed temperature, T = 0.10, and another at fixed field, B = 0.20. 
The magnetic phase diagram of the similar model for zero baryon density is investigated in 
Ref. [3] . The phase diagram at finite density is explored in Ref. [n| with the approximation 
f(u) ~ 1. We found that for T > 0.10, this approximation is no longer valid. 

Our main results demonstrate that for a given magnetic field and moderate tempera- 
ture, the most preferred nuclear phase in the SS holographic model is the multiquark-pion 
gradient (MQ-y</?) phase provided that the density is sufficiently large. We also study the 
temperature dependence of the baryon chemical potential, the free energy, and the linear 
response of the pion gradient of the mixed MQ-y^ phase and show that they inherit the 
temperature dependence mostly from the SS background. 

Extremely strong magnetic fields could have been produced in many situations. The 
Higgs mechanism in the cosmological electroweak phase transition could create enormous 
magnetic fields in the region between two different domains with different Higgs vacuum 



expectation values [18[ which could play vital role in the phase transitions of the nuclear 
soup at later times. At the hadron and heavy ion colliders, colliding energetic charged 
particles could produce exceptional strong magnetic field locally. The local magnetic fields 
produced at RHIC and LHC are estimated to be in the order of 10 14 ~ 15 Tesla 19]. On 
the astrophysical scale, certain types of neutron stars called the magnetars could produce 
magnetic fields as strong as 10 10 Tesla 20]. 

This article is organized as the following. In Section 2, the setup of the deconfined SS 
model with additional baryon vertex and string sources are discussed. Main results are 
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elaborated in Section 3. Section 4 concludes the article. 



II. HOLOGRAPHIC SETUP OF THE MAGNETIZED MULTIQUARK PHASE 



The setup we will use is the same as in Ref. [15|, the Sakai-Sugomoto model with ad- 
ditional baryon vertex and strings (baryon vertex is introduced in Ref. 2l|, |22[). Starting 
from a 10 dimensional type IIA string theory with one dimension compactified into a circle 
which we will label x A . Two stacks of D8-branes and D8-branes are then located at distance 
L from each other in the x 4 direction at the boundary. This separation will be fixed at 
the boundary and it will play the role of the fundamental scale of our holographic model. 
Open-string excitations with one end on the D8 and D8 will represent quarks with different 
chiralities. In the background where the D8 and D8 are parallel, excitations for each chirality 
are independent and there is a chiral symmetry in the background and at the boundary. For 
background with connecting D8 and D8, chiral symmetry is broken and there is a chiral 
condensate. When the energy of the connecting configuration is minimal and there is no 
extra sources, we define the corresponding boundary gauge matter to be in a vacuum phase. 

Since the partition function of the string theory in the bulk is conjectured to be equal to 
the partition function of the gauge theory on the boundary, the free energy of the bound- 
ary gauge matter is equivalent to the superstring action in the bulk (modulo a periodicity 



factor) 23[. We turn on non-normalizable modes of the gauge field a^,a^,a^ (defined in 
units of i?D4/2vra / ) in the D8-branes and identify them with the vector potential of the mag- 
netic field, B (defined in units of l/27ra'), the gradient of the chiral condensate, VV 9 ; an d 
the baryon chemical potential, fi, at the boundary respectively. These curious holographic 
correspondence between the branes' fields and the thermodynamical quantities of the gauge 
matter at the boundary allows us to study physics of the strongly coupled non-Abelian gauge 



matter at finite density in the presence of the external magnetic field. E 
be added using other components of the gauge field on the D8-branes 
not consider such cases here. 



ectric field can also 



16 



24| but we will 
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The background spacetime of the Sakai-Sugimoto model is in the form 

ds * = r^y 2 (f(u)dt 2 + + dx 4 2 ) + y 2 (u^i + 

F w = ^—e 4 , e* = g s (-^) , R 3 m = 7rg s NJ 3 s , 

V4 V R-D4 J 

where f(u) = l — u T /u 3 , = 167r 2 i? 3 D4 T 2 /9. V4 is the volume of the unit four-sphere and 
€4 represents the volume 4-form. l s and g s are the string length scale and the string coupling 
respectively. R is the compactified radius of the x 4 coordinate. This radius is different from 
the curvature Rda of the background in general. The dilaton field is denoted by (f) which 
will be eliminated by the function of u as stated above. 

The direction of the magnetic field is chosen so that the vector potential is 

4 = Bx 2 . (1) 

The baryon chemical potential \x of the corresponding gauge matter is identified with the 
non-normalizable mode of the DBI gauge field at the boundary by 

fx = Oq(u — > 00). (2) 

The five- dimensional Chern-Simon term of the D8-branes generates another axial part of 
the U(l), a^ 4 , by coupling it with B and Oq . In this way, the external magnetic field induces 
the axial current 2 a associated with the axial field af. The non-normalizable mode of this 
field at the boundary corresponds to the response of the chiral condensate to the magnetic 
field which we call the pion gradient, VV 9 - External field causes the condensate to form a 
domain wall which can be characterized by the gradient of the condensate with respect to 
the direction of the applied field. Therefore the pion gradient also acts as a source of the 
baryon density in our gauge matter. 

Additional sources of the baryon density and the baryon chemical potential can be added 



to the configuration in the form of the baryon vertex and strings 
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7^ u cV f(u c ) + n s {u c - u T ) , (3) 
J\f dfx source (4) 



where n s = k r /N c is the number of radial strings in the unit of 1/N C . Since the radial strings 
could merge with strings from other multiquark and generate a binding potential between 
the multiquarks, this number therefore represents the colour charges of the multiquark in the 
deconfmed phase. It is interesting to note that when there is only string source representing 
quark matter, the quark matter becomes thermodynamically unstable under density fluctu- 



ations 



However, adding baryon vertex together with the strings makes the multiquark 



10| . The multiquark phase is even more 



configuration stable under the density fluctuations 
thermodynamically preferred than the xS-QGP when the density is sufficiently large and 
the temperature is not too high. 

With this setup, then the DBI and the Chern-Simon actions of the D8-branes configura- 
tion can be calculated to be 



Scs 



JV r du u^ 2 J l + ^Jl + /(u)«) 2 - (a' v y + /WAf, (5) 

Juc V U 

3 poo 

—-JV / du (d 2 al ciq af' — d 2 a^ a^' 'af) , (6) 
2 J Uc 



where M = N R 2 Di / '(67r 2 (27ro/) 3 ) defines the brane tension. The factor 3/2 in the Chern- 
Simon action fixes the edge effect of the region with uniform magnetic field as explained in 
Ref. Q. 

We can write down the equations of motion with respect to each gauge field 



j A --Bfi + 3B a v, (7) 



a 2 



y/1 + f(u)( a^ - (a'v y + f(u)u^xf 

Vu 5 + B 2 u 2 atf , 3 . n tj a 

- d Baf (oo) + 3Baf . 



y/1 + /H«) 2 - K) 2 + f{u)uH'i 2 
d, ja are the corresponding density and current density at the boundary of the back- 
ground (u — y oo) given by 

. SS Pr „ 



j tM (x, u — > oo) 



(9) 



(d,f A ). (10) 



(3 



In terms of the gauge fields, they are 



d 



3a 



Vu 5 + 



B 2 u 2 a'V 



y/i + 7 R(jgg - K ) 2 + f{u)v^! 

Vu 5 + B 2 u 2 f(u)a[ A 



(12) 



^l + fiuKa'fy-ia'Vy + fiuWx'i 
In order to solve these equations, we need to specify the boundary conditions. Due to the 
holographic nature of the background spacetime, the boundary conditions correspond to 
physical requirement we impose to the gauge matter. If we want to address chirally broken 
phase of the gauge matter, we will take a^(oo) = \jip to be an order parameter of the 
chiral symmetry breaking and minimize the action with respect to it. This results in setting 
]a = 0. On the other hand, if we want to study the chiral-symmetric gauge matter (or 
chiral-symmetric quark-gluon plasma for N c = 3 and a^(oo) will be set to zero. 

For vacuum phase, Oq , a± and d,jA will be set to zero. 

In any cases, since the total action does not depend on x±{u) explicitly, the constant of 
motion gives 



u 3 [f(u)(C{u) + D(u) 2 ) - {j A - \Bn + 3Ba 



F 2 



where 



u 3 f(u 



\B{i + 3Ba%(u c )) 2 <K) 



(13) 



(14) 



'l + f(u c )u 3 c x'£(u c ) 

and C(n) = u 5 + B 2 u 2 , D{u) = d+3Bcii(u) — 3B\jip/2. The calculation of x'^Uc) is described 
in the Appendix as a result from the equilibrium and scale fixing condition 



x' 4 (u) du = 1. 



(15) 



The equations of motion Eqn. (Ej),® can be solved numerically under the constraint (TIB"]) . 
The value of n,\/<p,u c and the initial values of a^{u^), af (u c ) are chosen so that a^(oo) = 



fi, a^(oo) 



Vv 9 an d Lq = 1 are satisfied simultaneously. 



III. MAGNETIC PHASE DIAGRAM OF THE DENSE NUCLEAR PHASE 



Generically, the action ([5]) and are divergent from the u — > oo limit of the integration 
and we need to regulate it using the action of the vacuum which is also divergent. The 
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contribution from the region u — > oo is divergent even when the magnetic field is turned off 
and it is intrinsic to the DBI action in this background. The divergence can be understood 
as the infinite zero-point energy of the system and thus could be systematically removed by 
regularisation. 

Therefore the regulated free energy is given by 

J^ = n(fi 7 B)+fjui 7 (16) 

where Q,(/i,B) = S[a (u), a\{u)](e.o.m.) — S [magnetized vacuum]. Note that we need to 
Legendre transform the DBI and the Chern-Simon action to obtain the bulk action as a 
function of the non-normalizable modes a^af in order to identify it with the free energy 
of the gauge matter at the boundary. In terms of the free energy at the boundary, this is 
equivalent to the Legendre transform of the grand canonical with respect to fi and d. 

We can calculate the total action satisfying the equation of motion 
S[ao(u), ai(u)](e.o.m.) = Sds + Scs to be 

c _ K r r, r < \ I f(u)(l + f(u)u*x») 

Sm ~ N L { 1 f(u)(C(u) + Dluf) ~ (J a - IB, + SEW 

3 ^ ( a V( 3A - \Bn + 3Ba%) - f(u)D(u)af) J jfa + u^x'i 
Scs = -M-B / du ± Jyn> = ■ (18) 

sJf{u){C{u) + D{uf)- (j A -lBn + 3B a v) 

The three nuclear phases above the deconfinement temperature are governed by the 
same equations of motion, each with specific boundary conditions as the following, 

magnetized vacuum phase: Oq , af = 0; d, Ja = 0, 

multiquark-pion gradient phase: a^{u c ) = fi SO urce,ci^(u c ) = 0,a^(oo) = \/<fi,jA = 0, 
xS-QGP phase: x' 4 (u) = 0,a^(w c = ut) = 0, (oo) = 0,j A = f-B/i. 



We will demonstrate later that in the mixed phase, the pion gradient is generically dom- 
inated by the multiquark when the chiral symmetry is broken. In Ref. 15[, it is shown that 
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the pure pion gradient phase is always less preferred thermodynamically than the mixed 
phase of MQ-y^- It is interesting to note that for the pure pion gradient phase, a large 
magnetic field is required in order to stabilize the generated domain wall 13|. This critical 
field is determined by the mass of the pion in the condensate, B crit ~ m| / e. In Ref . 15[ , this 
critical behaviour is confirmed in the holographic SS model (the zero-temperature situation 



is studied in Ref. 



25[). More investigation of the pure pion gradient phase in the holographic 



model should be conducted especially when the field is large since the distinctive feature of 
physics from the DBI action becomes apparent in this limit. We will leave this task for 
future work and focus our attention to the mixed MQ-y^ phase in this article. 
The action of the magnetized vacuum when we set a^, = and d,j^ = is 



S [magnetized vacuum] 



U.(| 



yJC(u)(l + f(u)u*x2) 



du. 



where 



x' 4 (u)\ % 



f(u)u l 



3 ( /(u)u3C(u) 
f(u )u$C{u ) 



- 1 



(19) 



The position uq where x' A — > oo of the magnetized vacuum configuration increases slightly 
with temperature as is shown in Fig. [TJ The difference between each temperature decreases 
as the magnetic field gets larger and all curves converge to the same saturated value uq = 1.23 
in the large field limit. 

We can study the temperature dependence of the magnetized multiquark nuclear matter 
by considering its baryon chemical potential and the free energy as shown in Fig. [2j Both the 
chemical potential and the free energy decrease steadily as the temperature rises, regardless 

of the magnetic field. This is originated from the temperature dependence of f(u) = 1 J 

of the SS background in the deconfined phase. The temperature dependence could be fit 
very closely with the function 



(|r) 6 as the following 



/x = Li (d,B) x /l-(— )6, 



T 



F = F (d,B) 



T 
Jo 



(20) 
(21) 
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where for d = 1,B = 0.10; fio = 1.1849, Fq = 0.7976 respectively. For the baryon chemical 
potential (free energy), the best-fit value of To is 0.269 (0.233). The fittings are shown 
in Fig. [3J This could be explained by noting that the regulated free energy is given by 
fid + Q(fj,,B). The contribution from the first term is dominant therefore the free energy 
has almost the same temperature dependence as the chemical potential. However, there is a 

3 

minor contribution from Sbs + Scs containing f(u c ) = 1 — which for small temperature 
fractions modifies the temperature function in the following manner, 




n-«f-w (22) 

where are some arbitrary constants and Co, To are given by 

Co = d + C 2 , (23) 
■I = -J— (£ + ^1 (24) 

It should be noted from Fig. [3] that the temperature dependence is significant for T > 
0.10 and the approximation f(u) ~ 1 is not accurate for temperature in this range. The 
characteristic temperatures we found here are consistent with the phase diagram of the 
multiquark in Fig. [7J 

In the multiquark phase when the magnetic field is turned on, the pion gradient is induced 
by the field in addition to the multiquark. The multiquark phase thus contained the mixed 
content of multiquarks and the pion gradient. For moderate fields (not too large), the 
response is linear \/ip oc B. In contrast to the case of pure pion gradient phase, the domain 
wall in the mixed MQ-yv 9 phase is stable among the surrounding multiquarks even for small 
field. The critical magnetic field to stabilize the domain wall in the case of pure pion gradient 
is not required in the mixed phase. 

Figure[2](b) shows a linear relation between the pion gradient and the magnetic field which 
is valid up to moderate fields. For d = 1, we found that the slope, m (or the linear response), 
of this linear function depends on the temperature approximately as m = m Wl — (jr) 6 , 



and 




W ~ Bmoi l- - , (25) 
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where mo = 0.347, To = 0.177. The curve fitting is shown in Fig. HJ The density dependence 
is encoded in mo = mo(d),To = To(d). As the density increases, the slope of the linear 
response of the pion gradient becomes smaller as is shown in Fig. The ratio of the 
pion gradient density and the total baryonic density -R w = d^^/d = 3B\/(p/2d [l4] for 
B = 0.10, T = 0.10 is plotted in the log-scale in Fig. (b). It could be well approximated 
by 

R V(p ~ (const. )eT 6/5 , (26) 
3B 2 m 




2d »- > - < ■ < 27 > 
from Eqn. ( )25l) . This implies that the multiquark states are more preferred than the pion 
gradient in the presence of the magnetic field, the denser the nuclear matter, the more stable 
the multiquarks become and the lesser the population of the pion gradient. 

Finally we compare the free energy of the MQ- yy? phase and the chiral-symmetric quark- 
gluon plasma phase. For high density, d = 100, this is shown in Fig. [61 For a given density, 
the multiquark phase is more thermodynamically preferred than the xS-QGP for small and 
moderate fields. As the magnetic field gets larger, the xS-QGP becomes more thermody- 
namically preferred. When the field becomes very strong, the transition into the lowest 



Landau level finally occurs 26|]. For a fixed density, increasing magnetic field inevitably 
results in the chiral symmetry restoration. The phase transition between the MQ-y^ and 
the xS-QGP is a first order since the free energy is continuous at the transition and the 
slope has a discontinuity. It implies that the magnetization, M(d, B) = — , of the nuclear 
matter abruptly changes at the transition. 

On the other hand, for a fixed field and the moderate temperature, the increase in the 
baryon density could make the multiquark phase more stable than the xS-QGP. This is 
shown in the phase diagram in Fig. [7J At a given magnetic field, the multiquark phase 
could become the most preferred magnetized nuclear phase provided that the density is 
made sufficiently large and the temperature is not too high. In contrast, the effect of the 
temperature is the most dominant for chiral-symmetry restoration even when the field is 
turned on. Sufficiently large temperature will induce chiral-symmetry restoration for most 
densities as is shown the Fig. [7J(b). Nevertheless, theoretically we can always find sufficiently 
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large density above which the multiquark phase is more preferred. 

The transition line between the MQ-y and the xS-QGP phases in the (d, B) phase 
diagram can be approximated by a power-law 

B ~ d 0A38 (°- 436 ) (28) 

for the multiquark with n s = (0.2). This power-law is weaker than the transition line of the 
xS-QGP to the lowest Landau level studied in Ref. 26[ for the antipodal SS model (B ~ d 2 ^ 3 ). 
The multiquarks with more colour charges (n s > 0) are less preferred thermodynamically 
but they require higher densities. On the other hand, the transition line in the (d, T) 
phase diagram is an increasing function of d but weaker than the logarithmic of the density. 
Nevertheless, theoretically for a fixed B, T, we can always find sufficiently large density above 
which the MQ-yy? phase is preferred. The high density region is actually dominated by the 
multiquark phase indeed. 



IV. CONCLUSION 



We explore the properties of the miltiquark-domain wall (MQ-y<^) solution of the SS 
model above the deconfinement temperature. The temperature dependence of the baryon 
chemical potential, the pion gradient linear response (m), and the free energy of the MQ- 
SJf phase has been studied and fitted with a simple function, yfl — inherited from 
the deconfined SS background. Their characteristic temperatures, To, are different from 
one another depending on other parameters such as u c , the position of the baryon vertex. 
Remarkably, they do not depend on the field for moderate field strength B = 0.05 — 0.15. 

For chirally broken deconfined nuclear matter in the presence of the magnetic field, the 
nuclear matter with finite baryon density and chemical potential could respond to the mag- 
netic field by inducing a pion gradient or a domain wall of the chiral condensate. This pion 
gradient response is found to be a linear function of the field for moderate fields at any 
density. However, we demonstrate further that the population ratio of the pion gradient 
decreases as the density increases. The other sources of the baryon charge namely the mul- 
tiquarks finally dominate the chirally broken nuclear phase and most of the baryon density 
is in the form of the multiquark at high density. 
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Magnetic phase diagram of the dense gauge matter have been explored in the deconfmed 
SS model. At fixed magnetic field and moderate temperature, the MQ-yv 9 phase are more 
preferred than the xS-QGP for the high density region. The transition line in the (d, B) 
phase diagram at T = 0.10 can be fitted closely with the power-law B ~ d 0A38 (°- 436 ) for the 
multiquark with n s = (0.2). On the other hand, the transition line in the (d,T) phase 
diagram is weaker than the logarithmic of the density but nevertheless it is an increasing 
function with respect to the density. These imply that for sufficiently large density, the 
chirally broken multiquark phase is the most preferred nuclear phase even in the presence of 
the external magnetic field. 

The situation when density becomes extremely large and being dominant occurs in the 
core of dense star such as the neutron star. Therefore it is very likely that the core of dense 
warm star composes primarily of the multiquark nuclear matter even when an enormous 
magnetic field is present such as in the core of the magnetars. It is possible that a large 
population of the warm magnetars has multiquark cores. These warm dense objects could 
be relatively more massive than typical neutron stars. 
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Appendix A: Force condition of the multiquark configuration 

Fixing the characteristic scale Lq to 1 for the brane configuration requires balancing three 
forces in the gravity picture. The D8-brane tension must be in equilibrium with the tidal 
weight of the D4 source and the string tension of the colour strings. The derivation of the 
x' 4 (u c ) presented here is the same as in Ref. [ijj], it is included for completeness. 

We vary the total action with respect to u c to obtain the surface term. Imposing the 
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scale- fixing condition 2 dux' 4 (u) = L = 1, we found that 111 ] 



L(llr 



OS 



dx' A 



(Al) 



as the condition on u c . 

We perform the Legendre transformed action with respect to Oq' and a^' to obtain 



S 



L(x' 4 (u), d) du, 



Ml du , 

lu c V /(«) 



x Jf(u)(C(u) + D(u) 2 )- (j A --B/x + 35a 



n 2 

o ) ' 



(A2) 



where C(u) 



u ' 



B l u%B{u) = d + 35af(n) - 3Bsj(p/2. Note that the Chern-Simon 



action are included in the total action during the transformations. 
The Chern-Simon term with the derivatives a v ' a Al eliminated is 



.5', 



cs 



3 ,oc ( a V( jA - \Bn + Wa v ) - f(u)D(u)af) Ji+u*xf 
-N-B / du ± ;V ii_ . (A3) 



/(m)(C(m) + D{uf) -lj A - %Bu + Wa\ 



From Eqn. flAlj) .f lA2|) . (IA3j) .p|). and the boundary conditions, a^{u c ) = u source ,af(u c ) 
0, we can solve to obtain the condition for the static equilibrium 

1 r 9 (f c (C c + Dl) - (j A - \Bu + 3Ba%(u c ))^ 

{X A {U c )j 



fcUlVd? (1 + 1(^)3 + 3^^)2 

where f c = f(u c ),C c = C(u c ),D c = D{u c ). 
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FIG. 1: The position uq as a function of B of the magnetized vacuum for T = 0.02 — 0.15. The 
upper lines have higher temperatures. 
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FIG. 2: The chemical potential (a), the pion gradient (b) and the free energy (c) of the multiquark 
phase with baryon density d = 1 as a function of B for temperature T = 0.02 — 0.15. The lower 
curves represent multiquark at higher temperatures. 




FIG. 3: For d = 1, B = 0.10, (a) the baryon chemical potential as a function of T, the best-fit curve 
is in the form [iQyJl — (^) 6 with /j>q = 1.1849, To = 0.269; (b) the free energy as a function of T, 
the best-fit curve is in the form F ^Jl - (J) 6 with F = 0.7976, T = 0.233. Other curves within 
the range B = 0.05 — 0.15 can also be fitted well with the same Tq. 
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FIG. 4: The linear response or slope of the linear function between the pion gradient and the 
magnetic field as a function of the temperature for the range B = 0.05 — 0.15 and density d = 1. 



The red line is the best-fit curve in the form mo\/l — (^) 6 with mo = 0.347, To = 0.177 
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FIG. 5: (a) The pion gradient as a function of B for density d = 1, 10, 100 at T = 0.10. (b) The 
density ratio of the pion gradient with respect to the total baryon density of the multiquark phase 
at B = 0.10, T = 0.10 in the double-log scale. 
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FIG. 6: For the dense multiquark with d = 100, T = 0.10, (a) the chemical potential, (b) the free 
energy as a function of B. The multiquark curves in red are compared with the xS-QGP curves in 
blue for the chemical potential and the free energy. 



B 




0.16 
0.14 
0.12 
0.10 
0.08 
0.06 



20 40 60 80 100 













Afsb (n s = 


0)______ 








= 0.2) 




2 


4 6 


8 10 



FIG. 7: The phase diagram of the dense nuclear phases involving multiquarks when gluons are 
deconfined for (a) T = 0.10 and (b) B = 0.20. The chiral-symmetric quark-gluon plasma and the 
chirallly broken MQ-yy? phase are represented by xS and xSB respectively, n s is the number of 
colour strings in fractions of 1/N C . 
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